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ABSTRACT 

Modeling the structure formation in the universe, we extend the spherical 
collapse model in the context of MOND starting with the linear Newtonian struc- 
ture formation followed by the MONDian evolution. In MOND the formation of 
structures speed up without a need for dark matter. Starting with the top-hat 
over-dense distribution of the matter, the structures virialize with a power-law 
profile of the distribution of matter. We show that the virialization process takes 
place gradually from the center of the structure to the outer layers. In this sce- 
nario the smaller structures enter to the MONDian regime earlier and evolve 
faster, hence they are older than larger structures. We also show that the viri- 
alization of the structures occur in the MONDian regime, in which the smaller 
structures have stronger gravitational acceleration than the larger ones. This 
feature of the dynamical behavior of the structures is in agreement with this fact 
that the smaller structures as the globular clusters or galactic bulges have been 
formed earlier and need less dark matter in CDM scenario. 

Subject headings: gravitation-galaxies: formation-cosmology: theory-dark matter- 
large scale structure of universe 
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Introduction 



The conjecture for the existence of dark matter dates back to Zwicky in 1933, who failed 
to explain the dynamics of C oma cluster by the virial theorem through the distribution of 
visible matter (jZwicky 19331 ). To interpret the dynamics of the structure, the concept of 
missing mass or dark matter entered the Astrophysical studies since then. In addition to 
the cosmological scales, in the smaller galactic scales, t he flat rotati on curve of the spiral 
galaxies requires the existence of the dark matter halo (IBosma 19811 ). On the other hand 
studying the dynamics of the unive rse as a whole revea ls that the universe is dominated by 
the dark matter and dark energy (jSpergel et al. 20031 ). This hypothetical matter neither 
emits light nor interacts with the ordinary matter and only shows its presence through its 
gravitational interaction. The advantage of the Cold Dark Matter [CDM) model is that, 
it can successfully explain the rotation curve of the spiral galaxies and lensing by galaxies 
and cluster of galaxies. Within the framework of general relativity (GR) it also provides 
a reasonable description for the hierarchy in the structure formation. Although currentl y 
CDM and AC DM models are remarkably success ful in large scales (ISpergel et al. 20031 ). 
they cannot explain TuUy-Fisher and Freeman laws (IBosch fc Dalcanton 20001 ). On the other 
hand, the high-resolution N-body simulations are still in contradiction with the observations 
on sub-galac tic scales where th e y predict orders of magnitude more substructures than what 
is observed (IMoore et al. 19991: iKlvpin et al. 19991) . They also provide incompatible spatial 
distribution of the sub-halos (IKroupa et al. 20041 ). 



The other alternative to the dark matter is the modification to the gravity law by means 
of taking a generic form of action for the gravity rather than that of the Einstein-Hilbert 
action. This approa ch has been introduced to be an alternative model to the dark energy 
( ICarroU et al. 20041). f (R) gravity also is used to interpret the rotation curve of the spiral 
galaxies (ISobouti 20071 : ISaffari fc Rahvar 20081 ). There is also recent efforts on modifying 
the gravity law by using a simple kinetic Lagrangian whether the pressure can bend space- 
time sufficiently to replace the roles of dark ene rgy, cold da rk matter and heavy neutrinos in 
explaining anomalous accelerations at all scales (jZhao 20071 ). Halle et al (2008) also proposed 
a generalized lagrangian with non-uniform cosmological constant for the vacuum field within 
the framework of the Einstein gravity. 

Finally the third approach, which we are concerned with in this work is the modification 
of the conventional Newtonian law is so-called MOdified Newtonian Dynamic s (MOND). Th e 
dynamics of a structure in MOND under the gravitational field is given by (IMilgrom 19831 ): 



/^(^/«o)g = gN, (1) 

where qn is the Newtonian gravitational acceleration, = 1.2 x 10~^''ms~^ is a fundamental 



- 3 - 



acceleration parameter and fi{x) is a function for the transition from the Newtonian to the 
MONDian regime (e.g. fi{x) = x/vT+s^). The dynamics of the structure for g < 
deviates from the Newtonian law by and for <^ g we recover the Newtonian dynamics 
(i.e. /i = 1). On the other hand for g <^ oq, so-called deep MOND regime, ii{x) = x and the 
effective acceleration is given hj g = ^gN^o- 

Due to confusion in the definition of the dynamical concepts in MOND, this model can 
be interpreted as a modification to the gravity law instead of dynamics. Bekenestein and 
Milgrom (1984) used a non-conventional Newtonian action for the gravity to extract the 
modified Poisson equation as follows: 



V. [^(V0/ao)V0] = AnGp, 



(2) 



where in the spherical symmetric systems this equation reduces to equation ([T]). One of 
the problems with MOND is that it is not a covariant gravity model. Bekenestein (2004) 
proposed a covariant formulation of this model. This theory in addition to the metric has 
a scalar field as well as a 4-vector field called TeVeS, where in the limit of non-relativistic, 
small accelerations, the field equation reduces to that in MOND. The total potential in this 
theory is given by the sum of the Newtonian potential (pN and the potential due to the scalar 
field, (ps- 

$ = $jv + 0s, (3) 

where the added scalar field plays the role of the dark matter. TeVeS theory has also a 
Newtonian limit for the non-relativistic dynamics with significant acceleration. In the non- 
relativistic limit we can simplify the gravity inside the spherical system as 



9[r} 



9N, 



9n ^ oo; 
otherwise. 



(4) 



The advantage of MOND is that it could provide a successfu l fit to the rotation curve 



of spiral gala xies and dispersion velocity of elliptical galaxies (IMilgrom fc Sanderes 2003 



Sanders 19961 ). It has been tested against the Cosmic Background Radiation (jSkordis et al. 20061 ) 
gravitational lensing ( Chen and Zhao 20061: Zhao et al. 2006 : Angus et al. 2007h. stellar sys - 
tems and gal actic dynamics (Haghi et al. 20061; Nipoti et al 200?!: Tiret &: Combes 2007 ). 
solar system ( Bekenstein and Maeueiio 20061: Sereno and Jetzer 20061 ) and TuUy-Fisher and 
Freeman laws (IMcgaugh and Blok 19981 : iMcGaugh et al. 20001). MOND also decreases the 
mass discrepancy in the cluster of galaxies (jPointecouteau and Silk 20051 ) but yet in clusters 
it remain s necessary to invoke the undetected matter, possibly in the form of a mass ive 
neutrino (ISanders 19991 : ISanders 20031 : lAguirre et al. 20021 : ISanders and McGaugh 2002 I ) . 



Studying the cosmology and formation of the large structures in the universe is another 
tool to examine MOND. Studying the MONDian scenario of the structure formation has been 
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started by Felten (1984) and Sanders (1998). In the paper by Sanders (1998), it is shown that 
a patch of universe in the MONDian regime smaller than the horizon size can evolve with 
a different rate than the background, hence the structures naturally can be formed through 
this scale dependent dynamics. In this scenario small structures form before the larger ones 
and this provides a bottom-up hierarchical procedure for the formation of the structures in 
the universe. The problem with this model is that the center of collapse is not identified and 
every point in the space depending on the choice of the coordinate system can be considered 
as the center of the structure. For solving this problem one can consider MOND formula to 
be applied to a peculiar a c celeration dey eloping from density fluctuation rather than Hubble 



expansion (ISanders 20011 : iNusser 20021 ). Nusser (2002) used the amplitude of the CMB 
anisotropies as the initial condition in the N-body simulation to simulate the large scale 
structures in the universe. However to have a compatible result, one has to either decrease Qq 
by one order of magnitude or reduce the amplitude of the fluctuations at the initial condition. 
Recent studies in the formation of the galaxies by Sanders (2008) shows that massive elliptical 
galaxies may be formed at 2; > 10, as a consequence of the monolithic dissipation-less collapse. 
Applying MOND with cooling mechanism put an upper limit to the stellar clustering in the 
form of the galaxy. Extending the structure formation in TeVeS theory has been done by 
linear perturbation of metric, vecto r and scalar fields. The pr edictions are compatible with 



the observations of the structures (IDodelson fc Liguori 20061 ). Skordis (2008) also used a 



generalized TeVeS theory to construct the primordial adiabatic perturbations of a general 
family of the scalar field kinetic functions. 

In this work we extend the spherical collapse model in MOND for studying the general 
behavior of the structures during their formation. The initial baryonic density contrast for 
the structures is taken from the CMB anisotropy. We obtain the dynamics of the baryonic 
structures in the early epoches with the linear Newtonian structure formation until the 
entering of the structure to the MONDian regime and follow the evolution with the MOND. 
The dynamical evolution with MOND shows that the structures virialize with power-law 
distribution of matter. We show that while all of the structures re-collapse and virialize in 
the MONDian regime, the gravitational acceleration of the structures in this stage inversely 
depends on the size of the structure. This dynamical behavior of the structure formation 
is compatible with the less-existence of the dark matter in the globular clusters and central 
parts of the galaxies in CDM scenario. 

The organization of the paper is as follows: In section [2] we give a brief review of the 
spherical collapse model in MOND and extend it by looking at the dynamics of each layer in 
the onion model. In section [3] we discuss about the results of the calculation, showing that 
the density profile inside the structure is a time varying function during the evolution and 
the structure at the final stage, virialize from the center to the outer areas. In section H] we 
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summarize and discuss the results. 



2. Structure Formation in MOND 

In this section we introduce the MONDian cosmology and apply MOND in the spherical 
collapse model for studying the formation of the structures. Here we take the onion model 
for the spherical structure, dividing the sphere to the co-centric shells and studying the 
evolution of each layer separately and the structure as a whole. 



2.1. MONDian cosmology 

In the standard cosmology the dynamics of the universe in the matter dominated regime 
can be derived from the Newtonian gravity. Sanders (1998) used this approach to obtain 
the dynamics of universe in the MONDian scenario. In MODian cosmology, a patch of 
universe can evolve with a different rate than the background as soon as the acceleration 
fulfills the condition of g < oq. The result is dynamical decoupling of the smaller scales 
from the background which causes the production of the over dense regions in the universe. 
The reason for this feature of cosmological dynamics in MOND is that unlike the Newtonian 
mechanics, the acceleration in the comoving frame depends on the length scale. 

Let us take a spherical region with radius r from the background in which <C g{r). 
Using the Newtonian dynamics the acceleration is given by 

where M is the active gravitational mass and to have a compatible relation with the rel- 
ativistic results, we define it to be composed of the relativistic and non-relativistic matter 
as 

M=^(p + 3p), (6) 

where p and p are the density and pressure of the cosmic fluid. Substituting equation in 
equation ([5]), the acceleration is given by: 

r = -^^{p + 3p)r. (7) 

In equation ([7]), the gravitational acceleration increases linearly with r. This implies that 
there should be a critical radius Vc where inside it the acceleration is smaller than the 
MOND threshold oq and the dynamics is given by the MOND where outside that radius it 
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is Newtonian. This critical length scale obtain by equaling the left hand side of equation ([7]) 
with r = — as 

3ao /ON 



A7rG{p + 3p) 

This length scale separates the Newtonian (r > rc) and the MONDian (r < Vc) domains. As 
the density of universe changes with the expansion of the universe, the critical radius also 
changes with time. From the continuity equation, the matter and the radiation densities 
vary as p = poa~^ and p = poa~^. Using the definition of the critical density pc = 311^/8710, 
equation can be written in terms of the density parameters, Q, and the scale factor 

(9) 



A I 

where we adapt Ho = 75kms-^Mpc-\ = 0.02, ni°^ = 5 x lO'^ and Q^^ = 0. 

Now we do comparison of the size of a structure with r^. From equation the critical 
radius changes with the scale factor as oc in the radiation and oc in the matter 
dominant epoches. On the other hand the size of the structure is proportional to the scale 
factor (i.e. A oc a). So we expect that the Newtonian structures eventually will enter the 
MONDian regime as Vc grows faster than A. Using the adapted cosmological parameters, 
the critical radius at the present time is obtained — IOHq^ which means that the whole 
observable universe resides in a MONDian domain. Comparing these two length scales at 
the last scattering surface results in rc/H~^ ~ 10~^. The mass corresponding to this critical 
radius at last scattering surface is about Mc — and for M < the dynamics is given 

by MOND. While we expect to have density contrast growth for the scales r < Tc at the 
decoupling, comparing the critical mass Mc ~ WMq with the Jeans mass of Mj ~ IO^Mq 
at this time ind icates that the structures at these scales should be washed out by the pressure 



fISanders 1998[ ) 



The other feature of MONDian cosmology is that unlike the standard cosmology where 
decoupling redshift is smaller than the equality redshift {zde < Zeg), in MOND the equahty 
epoch is much after than the decoupling time. This feature results from this fact that 
decoupling is related to the baryonic density of the universe and the temperature and both 
parameters depend only on the scale factor, independent of the dynamics of the universe. 
So we expect to have the same decoupling redshift in the MOND as the standard cosmology. 
However since the dark matter does not exist in the MONDian cosmology, the equality will 
be shifted to the lower redshifts. For our adopted cosmological parameters the equality 
redshift is obtained Zeq = 400. 
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2.2. Structure Formation: Spherical Collapse Model 

In this part we model the evolution of a structure with an over-dense spherical region 
in MOND. To calculate the evolution of this spherical patch, for simplicity we take this 
over-dense region with a top-hat distribution of matter. As the acceleration of this structure 
depends on the distance from the center, we expect to have different dynamics for each 
radius. Hence, we divide the structure into the co-centric spherical shells like an onion 



model in cosmology (ILemaitre 19331 : iTolman 19341 : iBondi 19471 ) and calculate the dynamics 



of each shell separately. The initial density contrast for this over-dense region is taken from 
the fluctuations of the last scattering surface. In the Newtonian treatment of the structure 
formation (r^ < A), the density contrast grows linearly {S oc a) from the decoupling epoch 
up to the entrance of the structure to Vc- For the MONDian regime (A < Tc), we switch 
the dynamics to MOND an calculate the evolution of the structure. As an example let us 
take a sphere with a mass of M = lO^^M© and find its acceleration with respect to the 
center. In Figure ([T]) we compare the acceleration of this spherical structure in MONDian 
and Newtonian dynamics as a function of the redshift. We note that the redshift is defined 
according to the dynamics of the scale factor at the background. At the early epoches, the 
difference between these two dynamics is small as Vc < A, but after entering the structure 
to the MONDian regime A < r^, the evolution of the structure by the Newtonian dynamics 
and MOND start to diverge. This deviation of the dynamics from that of Newtonian plays 
the role of the dark matter in the standard scenario of the structure formation. 

For calculating the dynamics of each shell in the onion model, we take the following 
notation: r*(t) is the radius of ith shell as a function of time and tl^^ and rg„^ are the entering 
time and radius of the ith shell to the MONDian domain, respectively. The velocity of ith 
shell in terms of the Hubble parameter at entrance time is given by vl^^ = Hentf\nti^ ~ Knt)y 
where Hent is the Hubble parameter of the Newtonian background and ^enier is the density 
contrast of the sphere enveloped by r\ As we discussed in the previous section, all the shells 
will eventually enter the MONDian regime and we take this time as the initial condition for 
each shell in the MONDian evolution of the structures. Table ([T]) shows the initial density 
contrast, radius and the corresponding redshift of the entrance of each shell to the MONDian 
regime. In the MONDian regime, the acceleration approximately is y/gN^o and the evolution 
of each shell is as follows: 

(10) 



where M* is the mass of the structure enveloped by the ith shell. Using the initial conditions 
given by Table (1), we obtain the evolution of each shell as shown in Figure ([2]). To visualize 
the evolution of the shells, we divide the sphere into ten equidistant shells when all the 
structure is in the Newtonian regime and obtain their evolution as a function of background 
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redshift. The dynamics of shells shows that the inner shells evolve faster, reaching to a 
maximum radius and then collapse earlier than the outer ones. Here the initial radius of 
outermost shell is about lAkpc at the entrance time of Zenter ~ 146 to the MOND regime. 
This shell expands up to a maximum radius of ~ 49kpc at Zmax ~ 28. Eventually, the shell 
starts to collapse and virialize at z ~ 18 with the radius of ~ 28.5kpc. The maximum radius 
of each shell, Vmax, is obtained from integrating equation (fTOj) . letting r(t) = as follows: 

i 2 

The next phase of the evolution of shells after reaching to a maximum radius is re-collapsing. 
Similar to the standard scenario of the spherical collapse models we expect that the global 
radial velocity of the structure during the free fall collapse convert to the dispersion velocity 
and prevent the structure from a catastrophic collapse. This steady stage of the structure 
is given by the virial theorem. The corresponding radius that fulfill the virial condition is 
called the virial radius and is calculated from: 

-r'^^^ + V(r') =E, (12) 

where V{r^) is the gravitational potential at the ith shell. The potential in the Newtonian 
or MONDian regimes is given by 



^^"^^-^ OM^ + GAP+,/GW^,\n{rlJ. Newt. ^^^^ 

For the non-dissipative evolution of the structure, the total energy is conserved and we 
substitute the right hand side of equation (fT2l) by the energy of the system at the enterance 
time to the MONDian regime, E = ^H^ntKnt^ + ^/GM^\nrl^^. Using the potentials given 
by equation (fT3l) at the left hand side of the equation (fT2l) . the virial radius for the MOND 
and Newtonian regimes obtain as 

r^„,e"-i/2, MOND; 

rlnti^-^^^f^)-'- Newt. (14) 

ent 

If the virialization of the structure takes place in the MONDian regime, r*^^ > r*^^ implies 
a > 1/2 which results in 2T^^^/W^^^ > 1 where T^^^ is the kinetic energy for a unit mass 
and W^ent — 'r'^dV^/dr\ On the other hand if the structure virializes in the Newtonian 
regime, r*j^ < r^^^ condition from the equation f|T^ implies 2T^j^jW^j^^ < 1. We calculate 
the expression of '^T^nt/^ent esich shell (see Table [1]) and show that all the shells of the 
structure viriahze in the MONDian regime (i.e. r*j^ > rg^J. 
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Table ([2]) shows the parameters of shells for the moment of maximum radius and the 
virialization stage and Fig. ([2]) visualizes quantitative behavior of the shells during their 
evolution. In this figure the spots on the evolution curves shows the two critical stages of the 
maximum radius and the virialization radius, reported in Table ([2]). The evolution lines of 
the shells show that the inner shells evolve faster and virialize at the higher redshifts, while 
the outer shells evolve slower. 



3. Predictions of the Model 

In this section we discuss the evolution of the density contrast and the profile of matter 
distribution inside the structure. We use the definition of the density contrast of the shells 
in our model: 

5 it) = , (15) 

where is the density of ith shell and p is the density of the background. As the collapsing 
of the shells starts from the inner to the outer parts of the structure, we will not have a 
shell crossing during the evolution of the structure and from the conservation of the mass, 
we can perform the Jacobian transformation from the initial distribution of matter inside 
the sphere to the evolved distribution as follows: 

p' = pu^r^, (16) 

where represents the thickness of the ith shell. Substituting the dynamics of each shell 
from the previous section in equation (fT6ll . we obtain the evolution of density contrast for 
the shells up to the virialization stage as shown in Figure ([3]). The inner shells evolve faster 
and reach the non-linear regime [6 > 1) at the higher redshifts while the outer ones evolve 
slower. Figure (HI) shows the dependence of the corresponding non-linear redshift to the 
mass of structure. Comparing a small scale structure with the mass of ~ IO^Mq with a 
galaxy having the mass of ~ IO^^Mq shows that the former structure enters the non-linear 
regime at Zni — 106 while the later one becomes non-linear at Zni — 33. More details on the 
characteristic redshifts of the structures in terms of their masses is reported in Table ([2]). 
From this table we extrapolate the dependence of the non-linear redshift, maximum radius 
redshift and virialization redshift to the mass of structure with the following functions (see 
Fig. H: 

M 

\ogizni) = -0.1671og(— ) + 3.38, 

Mq 

M 

log(z„a.) = -0.1921og(— ) + 3.56, 

Mq 
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M 

log(z„,) = -0.2031og(— ) + 3.48. (17) 

Mq 

In what follows we describe the qualitative predictions of this simple model. 

3.1. Density Profile 

In this part we compare the evolution of the density profile of the structures in the 
spherical collapse model for the Newtonian and MONDian regimes. In the Newtonian regime 
{g > ao) we have seen that the dynamics of the shells in the structure depends only on time 
and transforming to a comoving frame, the dynamics is scale independent (see equation [7j) . 
This means that the dynamics is invariant under the scale transformation and the result 
is preserving the initial profile of the structure. In MOND {g < a^), since the dynamics 
depends on the scale unlike the Newtonian case the density profile will change with time. 
In Fig. ([5]) we plot the spatial variation of density from equation ( fT6|l for a galaxy mass 
structure in four different stages of z = 400, 87, 23 and 18. The initial stage at the redshift 
of z = 400 is taken a top-hat distribution for the density when the innermost shell enters to 
the MONDian regime. We continuously do Jacobian transformation from this stage to the 
later times (i.e. z < 400) and obtain the density of each shell. While the structure evolve, 
the distance of shells change with time as shown with a point on the profile representing 
the position of each shell (Fig. [5]). As the inner shells enter to the MONDian regime earlier 
we expect the density profile deviates from the homogeneous distribution. Fig. ([5]) shows a 
small deviation of the density profile from the homogenous one at z = 87. At this moment 
all the shells are in the MONDian domain, however as the inner shells have been evolved 
faster, they are more dense than the outer shells. Each shell after virialization freezes and 
preserves its density. We fit the density profile of the structure when the outermost shell, 
the latest part of the structure, virializes and the result is a power-law function of p oc 
with the index of /? ~ —1.22. 



3.2. Age of a structure 

A question that may be answered in this simple model is the dependence of the age of 
the structure to the mass. Observations show that small isolated structures as the glob ular 
clusters or the center of galaxies are older than the spiral arms (jChanoyer et al. 19981 ). In 
this model we have seen that the smaller structures enter the MONDian regime earlier, evolve 
faster and virialize at the higher redshifts compared to the larger structures. As an example. 



- 11 - 



in Table ([2]) it is shown that a structure with the mass of IO^Mq viriahzes at 2; ~ 71 while 
the whole of galaxy virializes at z ^ 18. The age of a structure from the virialization up 
to now is given by tane = C° dt, where to is the present age of the Universe. Changing the 
variable to the redshift results in: 



(1 + z)^d^\i + zf + ni'\i + zY + (1 - fif^)(i + zY 

where we have adopted the cosmological parameters in (12. ip . The age of structures depending 
on their masses is given in Table ([2]). 




3.3. Contribution of dark matter in the smaller structures 

The other question in the formation of the structures is that why the smaller structures 
are mainly made of the baryonic matter and have a less dark matter compared to the larger 
ones. This feature of the structures can be explained if we can show that in addition to 
the fast evolution of the smaller structures, the acceleration of the structure in terms of oq, 
{.Qn/clo) gets larger. To show this property of the structures in MOND, we calculate the 
gravitational acceleration at the virialization time of each shell and compare it with qq. As 
the density of a structure at this stage changes with p oc r^, the gravitational acceleration 
will depend on r as QNij) oc r^^^. For [3 = —1.22 the smaller radii should have larger 
acceleration compare to the larger ones. 

In Figure ([6]) we plot the gravitational acceleration for each shell at the time of viriliza- 
tion in terms of the virialization radius. For the outer shells the gravitational acceleration, 
QN/do is smaller than the inner shells. This means that the smaller structures after viri- 
alization need less dark matter compare to the larger ones in the standard CDM scenario. 
Finally we plot the mass of the structure after virialization in terms of its size in Figure 
([7]). This feature reveals the general property of structures in MOND and an accurate result 
to compare with the observation may be obtained by N-body simulation of the structure 
formation. 



4. Conclusion 

Summarizing this work, we extend the spherical collapse model to study the generic 
properties of the structure formation in MOND. We showed that in the MONDian scenario, 
structures can evolve without a need to the dark matter and have the following three main 
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features: (a) MONDian spherical collapse unlike to that of CDM does not preserve the initial 
density profile of the structure. In MOND starting with an initial homogenous density profile 
the structure virializes with a power-law distribution of the matter, having singularity at the 
center, (b) We showed that the small scale structures enter the MONDian regime earlier 
than the larger ones, evolve faster and virialize at higher redshifts. This picture from the 
spherical collapse model in MOND provides a bottom-top scenario for the structure formation 
in which the smaller structures formed before the larger ones. This result is compatible with 
the observations where the old stars are located in the smaller structures such as the globular 
clusters or center of galaxies, (c) Finally we showed that while the smaller structures enter 
the MONDian regime before the larger ones, they have larger gravitational acceleration at 
the virialization time and hence need a less dark matter in CDM scenario. 

We would like to thank anonymous referee for useful comments. Also we would like to 
thank M. Nouri-Zonoz for reading the manuscript and giving useful comments. 
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Fig. 1. — Comparison of the accelerations in the Newtonian (dashed-hne) and MONDian 
(sohd-hne) regimes for a region with a galaxy mass scale (M = IO^^Mq) as a function 
of redshift. For the early universe the difference between the two dynamics is small but 
increases at the later times. The corresponding redshift for entering of the structure to the 
critical radius is Zenter = 146. 



- 16 - 



50 
45 
40 
35 
30 

'o' 

^25 
20 
15 
10 

5 



2.1 2 1.9 1.8 1.7 1.6 1.5 1.4 1.3 1.2 

log (Z) 

Fig. 2. — The evolution of shells inside the spherical structure with a galaxy mass (lO^^M©) 
as a function of redshift, = ri{z). We take ten equidistant shells at the initial condition 
for representing the evolution of each shell. The trend of the evolution is reaching to a 
maximum radius, recoUapsing and virializing . The dashed line connects the maximum radii 
of the shells and the dotted line indicates the virialized line, connecting the virilize points of 
each shell. 
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Fig. 3. — A log-log plot of the density contrast evolution for each shell as a function of the 
background redshift. Curves from the left to the right corresponds the inner to the outer 
shells of the structure. The horizon line represents S = 1, separate the linear and non-linear 
regimes. The inner layers reach to the non-linear regime at the higher redshifts while the 
outer ones become non-linear at the lower redshifts. 
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Fig. 4. — Dependence of the non-linear redshift (solid line), maximum radius redshift 
(dashed-line) and virialized redshift (dotted-dashed hne) as a function of the mass of struc- 
ture in logarithmic scale. Smaller structures enter to the non-linear regime at higher redshifts 
and the larger ones enter at lower redshifts. This feature of dependence of mass to the char- 
acteristic redshifts reveals that the smaller structure in this model form before the larger 
ones. 
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Fig. 5. — The evolution of density profile in logarithmic scale normalized to the initial 
density y = p/pmiuai, as a function of distance from the center of structure in four different 
redshifts. The initial density profile is taken top-hat distribution with ten equidistant shells 
for representing the evolution of the density. The inner most shell enters to the MOND 
domain at z = 400. Each point on the curves notifies the position of the shell. We plot 
density profile of the structure for z = 400 and z = 87 aX the left panel and z = 23 and 
z = 18 at the right panel for comparison. At z = 17 all the shells virialize and the density 
profile freezes with a power-low function of p oc where (3 ~ —1.22. 



- 20 - 




Fig. 6. — The Newtonian gravitational acceleration of each shell normalized to as a 
function of shell size at the virialization time. Spots on the curve represent the acceleration 
and the position of each shell. The gravitational acceleration in the inner shells is larger 
than the outer shells. 
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Fig. 7. — The mass of structure in terms of size of structure after virialization. Spots on the 
curve represent the the corresponding value for each shell. Structures with the mass smaller 
than ~ IO^^Mq are virializad with the radius smaller than ~ lOkpc. 
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Table 1: The parameters of each shell at the entrance time to the MONDian regime. The 
first column represents the index of shell. The shells are taken equidistant when all the 
structure is in the Newtonian regime. The second column indicates the mass that enveloped 
by the corresponding shell. The third column shows the radius of each shell at the time of 
crossing the MOND domain and the forth column is the corresponding redshift. The fifth 
column is the density contrast at the entrance time and the sixth column shows the ratio of 
the kinetic energy per mass to W — rdV/ dr of each shell at the entrance time. 

~l M^lO^^o] n{Unter)[kpc] Z.neter Sinter X l^''' ^^iJ^l,^ 



1 0.001 0.45 400 2.71 2.24 

2 0.008 1.25 297 3.70 2.28 

3 0.027 2.40 242 4.52 2.32 

4 0.064 3.55 218 5.03 2.36 

5 0.124 5.03 195 5.61 2.41 

6 0.215 6.35 184 5.97 2.48 

7 0.341 8.17 170 6.47 2.51 

8 0.510 9.77 161 6.81 2.53 

9 0.725 11.44 154 7.12 2.55 

10 1.000 13.57 146 7.51 2.57 
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Table 2: The parameters of shells in the spherical collapse model at the maximum radius and 
viriahzation redshift. The first column represents the index of shells. The second column 
indicates the mass that enveloped by the corresponding shell. The third, forth and the 
fifth columns correspond to the maximum radius, redshift and the density contrast at that 
moment, respectively. The sixth, seventh and the eighth columns show the virialize radius, 
redshift and the density contrast, respectively. The ninth column indicates the redshift that 
a shell enters to the non-linear regime {6 > 1) and the last column is the age of the structure 
after viriahzation. 
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0.001 


1.32 


106.52 


1.01 


0.79 


71.46 


28.51 


106.52 


12.48 
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0.008 


3.72 


73.94 


1.25 


2.21 


49.84 


34.37 


78.57 


12.43 
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0.027 


7.43 


58.88 


1.44 


4.45 


38.52 


36.82 


65.66 


12.38 
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0.064 


11.24 


49.53 


1.52 


6.85 


32.57 


38.32 


54.28 


12.34 
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0.124 


16.24 


43.44 


1.65 


9.79 


28.41 


40.30 


48.75 


12.29 


6 


0.215 


20.77 


39.26 


1.72 


12.46 


25.57 


42.55 


44.58 


12.25 
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0.341 


27.72 


35.23 


1.84 


16.01 


22.80 


44.49 


40.66 


12.19 
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32.65 


1.89 
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21.11 


45.05 


38.19 


12.17 
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0.725 


38.63 


30.75 


1.96 


24.01 
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36.32 


12.15 


10 


1.000 


16.98 


28.02 


2.06 


28.51 


17.99 


18.72 


33.66 


12.09 



